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ABSTRACT 


The  effect  of  the  central  obscuration  of  an  annular  pupil  of  an  incoherent 
optical  imaging  system  on  its  aberrated  point-spread  function  (PSF)  is  discussed 
A  simple  model  is  presented  which  approximates  the  aberrated  PSF  by  an 
aberration-free  PSF  scaled  by  its  Strehl  ratio.  The  approximation  of  a  PSF  with 
rotationally  symmetric  aberrations  by  an  aberration-free  PSF  scaled  by  the  Strehl 
ratio  improves  as  the  obscuration  of  the  pupil  increases.  Specifically,  as  the 
obscuration  approaches  the  outer  radius,  the  PSF  of  optical  systems  with 
rotationally  symmetric  aberrations  becomes  increasingly  similar  to  the  aberration- 
free  PSF  scaled  by  the  Strehl  ratio  in  regions  within  the  central  Airy  disk  A  PSF 
with  balanced  coma  aberration  can  also  be  approximated  with  similar  accuracy 
and  over  the  same  region  as  for  rotationally  symmetric  aberrations  For  the 
rotationally  nonsymmetric  aberrations  of  astigmatism,  balanced  astigmatism,  and 
coma  the  aberrated  PSF  in  its  central  region  can  also  be  approximated  by  the 
scaled  aberration-free  PSF  but  for  large  Strehl  ratios.  For  rotationally 
nonsymmetric  aberrations,  the  region  or  spot  sizes  are  given  for  which  the 
aberration-free  encircled  energy  without  scaling  by  the  Strehl  ratio  provides 
better  agreement  with  the  actual  aberrated  encircled  energy  than  the  aberration- 
free  encircled  energy  scaled  by  the  Strehl  ratio. 


CHAPTER  1 


INTRODUCTION 

This  thesis  describes  a  simple  model  that  can  be  used  to  approximately 
calculate  the  point-spread  function  (PSF),  Strehi  ratio,  and  encircled  energy  for  an 
aberrated  imaging  or  a  laser  transmitter  system  This  work  is  an  extension  of  the  work 
done  by  Mahajan,1  who  developed  a  simple  model  for  estimating  the  effects  of 
rotationally  symmetric  aberrations  on  the  PSF  of  a  system  with  a  circular  pupil.  It  has 
been  shown1-2  for  circular  pupils  that  the  irradiance  distribution  within  the  central 
region  of  the  Airy  disk  of  a  rotationally  symmetric  aberrated  PSF  is  largely  unchanged 
from  the  aberration-free  case  except  for  a  scaling  factor  equal  to  the  Strehi  ratio.  The 
current  work  focuses  on  the  applicability  of  such  a  model  to  syst  ems  with  annular 
pupils,  for  both  rotationally  symmetric  and  nonsymmetric  aberrations,  with  uniform 
and  Gaussian  illuminations.  Gaussian  illumination  is  of  interest  for  example,  to 
determine  the  focal-plane  distribution  of  a  TEM0  0  mode  laser  beam. 

The  irradiance  distribution,  or  PSF,  and  the  encircled  energy  for  an  aberrated 
system  with  an  annular  pupil  can  also  be  estimated  from  a  knowledge  of  the  Strehi 
ratio  and  the  aberration-free  PSF.  This  holds  true  for  defocus,  spherical,  balanced 
spherical,  and  balanced  coma  aberration.  They  can  also  be  estimated  for  astigmatism, 
balanced  astigmatism,  and  coma  but  for  smaller  spot  radii  and  larger  Strehi  ratios.  In 
particular,  the  PSF  and  encircled  energy  for  an  aberrated  system  can  be  estimated  by 
scaling  the  aberration-free  PSF  and  encircled  energy  by  a  factor  equal  to  the  Strehi 
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latio  as  was  shown  for  circular  pupils  by  Ma’najan  The  difference  between  the  actual 
aberrated  PSF  (and  encircled  energy)  and  the  aberration-free  PSF  (and  encircled 
energy)  scaled  by  the  Streh!  ratio  decreases  as  the  obscuration  ratio  increases  for 
rotationally  symmetric  aberrations  The  percentage  difference  between  estimated  and 
actual  encircled  energy  is  less  than  20%  within  the  Airy  disk  for  Strehl  ratios  S  >  0  4 
for  rotationally  symmetric  aberrations,  balanced  coma,  and  astigmatism  for  all 
obscurations  The  percentage  difference  drops  to  less  than  5%  for  rotationally 
symmetric  aberrations  for  Strehl  ratios  as  low  as  0  1  for  an  obscuration  ratio  of  0.9. 

The  difference  between  the  actual  and  estimated  encircled  energies,  in.  general, 
increases  (up  to  an  asymptotic  level)  as  the  spot  radius  (used  to  calculate  encircled 
energy)  increases.  This  is  a  result  of  encircled  energies  of  aberrated  PSFs  converging 
to  the  unaberrated  values  as  the  radius  increases.  The  region  or  spot  sizes  for  which 
the  aberration-free  encircled  energy  without  scaling  by  the  Strehl  ratio  provides  better 
agreement  with  the  actual  aberrated  encircled  energy  than  the  aberration-free  encircled 
energy  scaled  by  the  Strehl  ratio  depends  on  the  type  and  amount  of  aberration  The 
smallest  region  is  for  the  case  of  coma  and  is  approximately  half  the  Airy  disk  Other 
types  of  aberrations  ;  spot  sizes  ranging  from  the  Airy  disk  (for  rotationally 
nonsymmetric  aberrations  and  rotationally  symmetric  aberrations  with  unobscured 
pupils)  up  to  several  times  the  Airy  disk  (for  rotationally  symmetric  aberrations  with 
large  obscurations) 
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CHAPTER  2 


THEORY  AND  BACKGROUND 

The  irradiance  distribution  of  the  image  of  an  incoherent  point  object  of 
wavelength  A.  formed  by  a  system  with  an  aperture  or  pupil  is  called  its  diffraction 
point-spread  function  (PSF)  It  is  given  by3 

.«■  a  ^  2 

f  fG(^Tl)exP[-“(x^  +  >;Tl)ld^ri  (l) 

(M)  Li 

G(4.3l)  's  the  pupil  function  of  the  system  given  by 

G(£,,r|)  =  .4(£,rt)  exp[(2rti/X)ff(^,ri)]  at  points  inside  the  pupil  (2) 

=  0  at  points  outside  the  pupil, 

where  A(^,r ))  is  the  amplitude  and  W^.r))  is  the  wave  aberration  at  a  point  (£,r|)  on 
the  pupil  W^ti)  is  the  optical  path  length  difference  between  the  aberrated 
wavefront  and  the  Gaussian  reference  sphere.  The  reference  sphere  passes  through  the 
center  of  the  exit  pupil  and  has  a  radius  of  curvature  R.  It  is  centered  at  the  Gaussian 
image  point,  the  origin  of  the  (x,y)  image  plane 

The  diffraction  PSF  written  in  terms  of  polar  coordinates  and  located  in  a  plane 
normal  to  the  z  axis  at  a  distance  z  from  the  exit  pupil  (such  that  z  does  not  necessarily 
equal  R)  is,4-3 

PS  1 2n  r?  2 

/(r,G,;r;e)  =  --^-T  f  J ^(p)exp[/d>(p,0)]exp[— 7U — rpcos(0 - ©, )]pdpd0  .  (3) 

Here,  ^  =  pcosB,  rj  =  psinG,  x  -  rcosQ/,  y  =  rsinGy,  e  <  p  £  1  and  0  <  0  <  In.  <D(p,  0) 
is  the  phase  aberration  and  is  related  to  the  wave  aberration  according  to  <t>  = 
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(2n/X)lV.  The  explicit  dependence  on  z  is  of  interest  in  the  case  of  a  defocused  system 
where  z*  R,  otherwise,  z  =  /?,  and  Eq.  (3)  reduces  to  Eq  (1)  with  a  change  of 
coordinates  The  integral  is  taken  over  the  clear  region  of  the  pupil  with  coordinates 
(p,  0),  where  p  is  in  units  of  the  outer  radius  a  of  the  pupil  The  obscuration  ratio,  e, 
is  defined  such  that  at  is  the  inner  radius  of  the  pupil.  The  coordinates  of  the 
observation  point  of  the  PSF  are  (r,  0,,  z)  where  r  is  in  units  of  X  R  2a  -  X  F  (/•'  being 
the  f-number  or  the  focal  ratio  of  the  Gaussian  image  forming  light  cone).  The 
amplitude  A(p)  is  assumed  to  be  rotationally  symmetric  in  view  of  Gaussian  beams 
discussed  later.  The  total  power  transmitted  through  the  exit  pupil  is  P,  and  Sp  = 
tw:(1  -  e:)  is  the  area  of  the  exit  pupil. 

If  the  wave  aberration  is  expanded  into  a  power  series,  the  resulting  aberration 
terms  of  lowest  order  are  of  degree  4  and  are  called  the  primary  (or  Seidel) 
aberrations.6  The  primary  aberration  {unction  may  be  written  in  polar  coordinates 
as6'7 

W(p,  0)  =  As  p *+Ac  p3cos0+/4fl  p2cos20+/lar  p 2+A,  pcos0,  (4) 

where  A,  is  the  peak  value  of  the  corresponding  aberration  term  Thus  As  is  the  peak 
value  for  spherical,  Ac  is  for  coma,  Aa  is  for  astigmatism,  Aj  is  for  defocus,  and  A,  is 
for  tilt. 

The  aberration  function  may  also  be  expanded  in  terms  of  a  complete  set  of 
Zernike  circle  polynomials  in  the  case  of  a  circular  pupil*  and  Zernike  annular 
polynomials  in  the  case  of  an  annular  pupil  9  The  Zernike  polynomials  are  useful  in 
that  the  aberration  terms  are  of  the  form1011 

=  cm  y]2(n+  1)  /C(p:e)cosOT0,  (5) 

where  n  and  m  are  positive  integers  (including  zero),  n  -  m  >  0  and  even, 
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(6) 


em  =  1  /  V2 ,  m  =  0 
£ra  =  l,  /»*0, 

and  cnm  represents  the  standard  deviation  of  the  aberration  across  the  pupil  (unless  n  = 
m  =  0).  Each  Zernike  polynomial  is  composed  of  terms  of  the  form  A  pp  cos’  0  , 

where  p  and  q  are  positive  integers,  such  that  the  variance  of  the  aberration  is 
minimized.10  That  is,  an  aberration  of  a  certain  order  in  the  power  series  expansion  is 
mixed  or  balanced  with  aberrations  of  lower  order  such  that  its  variance  is  minimized 
thereby  maximizing  the  central  irradiance  for  small  aberrations.9- I0- 11  The  standard 
deviation,  0X4,;  e),  is  given  by12 

cj  =  (a>2)-(d>)2  (7) 

such  that  a,*,  =  cnm  if  <J>  is  expressed  as  an  orthonormal  Zernike  polynomial 
The  average  values  of  the  wth  power  of  are  defined  according  to 

1  2n  /  1  2* 

(<*>")  =  J  J^i(p)a>',(p,e)pdpde  /  j  J.4(P)pdPde.  (8) 

e  o  /  e  o 

The  primary  aberrations  and  their  standard  deviations,  o^,;  e),  have  been 
summarized  by  Mahajan13  and  are  shown  in  Table  1  for  the  case  of  uniform 
illumination  Note  that  Szapiel  has  studied  aberration  balancing  by  extending  the 
approach  used  by  Marechal  to  non-uniform,  radially  symmetric  amplitude 
distributions.14  He  finds  that  for  the  case  of  a  Gaussian  amplitude  (or  apodization) 
function,  Zernike  polynomials  result  in  optimal  aberration  balancing  when  the  Gaussian 
beam  waist  is  three  or  more  times  the  pupil  radius.  For  weakly  truncated  pupils  where 
the  pupil  is  about  three  or  more  times  the  Gaussian  beam  waist,  however,  the  optimum 
aberration  balancing  is  obtained  using  Laguerre  polynomials.14  Mahajan  has  followed 
the  same  Gram-Schmidt  orthogonalization  process  used  to  derive  Zernike  annular 
polynomials  to  generate  orthogonal  aberration  polynomials  for  Gaussian  illumination 
of  annular  pupils 9  Standard  deviations  of  aberration  functions  and  representations  of 
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Table  1 

Primary  Aberrations  and  Their  Standard  Deviations  for  Optica  1  Systems  with  Uniformly  Illuminated  Annular  Pupils 


balanced  aberrations  for  the  case  of  Gaussian  illuminated  annular  pupils  have  also  been 
reported  by  Mahajan.11 

The  ratio  of  central  irradiances  (i.  e.,  at  the  origin)  of  aberrated  and 
unaberrated  PSFs  is  called  the  Strehl  ratio  and  is  given  by13-15 

5  =  /( 0)*  /(0)*.0  (9) 

j  j  /I(p)exp[/O(p,0)]pdpd0 

C  0 _ 

j  j/j(p)pdpde 

E  0 

Marechal,16  Born  and  Wolf,17  Mahajan,13  and  Szapiel,14  among  others,  have  discussed 
approximate  expressions  for  the  Strehl  ratio  based  on  the  aberration  variance: 


Sm  =  {\-^V2Y 

(11) 

Sg  =  exp[-(o^], 

(12) 

where  the  subscripts  m  and  g  refer  to  the  Marechal  and  Gaussian  approximations 
respectively  Mahajan  has  shown  that  the  Gaussian  approximation  gives  a  better 
approximation  than  the  Marechal  approximation  for  classical  as  well  as  balanced 
(Zernike)  aberrations  in  systems  with  uniformly  illuminated  pupils.1318  Szapiel 
presents  a  method  for  modifying  the  Marechal  approximation  to  account  for  non 
uniform,  radially  symmetric  apodization  functions.  The  method  described  by  Szapiel 
does  not  require  an  explicit  analytical  expression  for  optimum  balanced  wavefronts, 
aberration  variance,  or  minimization  derivatives.14  Instead,  the  approach  requires  the 
calculation  of  moments  of  the  apodizing  function  and  resulting  determinants  of  a  set  of 
linear  equations  of  aberration  functions  which  minimize  the  aberration  variance 

Conversely,  the  maximum  aberration  for  a  given  Strehl  ratio  can  be  estimated 
and  is  referred  to  as  the  optical  tolerance  The  Marechal  tolerance  condition  for  a 


i 
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Strehl  ratio  of  greater  than  or  equal  to  0.8  at  the  diffraction  focus  is  that  the  root-mean 
square  departure  of  the  wave-front  from  the  reference  sphere  that  is  centered  on  the 
diffraction  focus  shall  not  exceed  the  value  XI 14;  i.  e.,  oM,  <  X  /14  which  follows  from 
S  =  0.8  >  l-(2;t  /  A,)2(aM,2).19  Although  not  explicit  in  the  Strehl  approximations 
above,  the  Marechal  tolerance  condition  is  relatively  insensitive  to  variation  of 
obscuration  ratio  and  the  type  of  aberration.  The  Rayleigh  criterion  on  the  other  hand 
specifies  that  for  a  Strehl  ratio  of  .S’  =  0.8,  fVrnax  =  X/4,  which  is  applicable  only  to 
primary  spherical  aberration20  for  e  s  0. 

The  PSF  for  rotationally  symmetric  aberrations  may  further  be  simplified  using 
the  identity21 

j  2r 

J0(x)~  —  Jexp[/xcos(0-a)]d0,  (13) 

and  so  the  PSF  in  Eq.  (3)  becomes  proportional  to  the  modulus  square  of  the  Hankie 

transform  of  the  pupil  function, 

4  PS  !>  d 

7(/-;  r;e)  =  — — ^  f >J(p)exp[/cI>(p)]y0(7r/'p — )pdp  (14) 

(M  i  - 

The  corresponding  fraction  of  the  total  energy  contained  within  a  circle  of  radius  rc 
(the  encircled  energy)  is 

rc 

E(rc)=  2njl{r)r  dr .  (15) 

0 


For  the  aberration-free  case,  <P(p)  =  0,  with  uniform  illumination,  /l(p)  =  1,  the 


PSF  in  the  plane  of  the  geometric  focus,  z  =  R,  becomes22 


I(r,R,z)  = 


1  PSp  \2Jfnr)  z2Jfnzr) 


(XRy  (I  -e‘)‘  I  nr 


In  the  case  of  a  circular  pupil,  e  =  0,  and  Eq.  (16)  reduces  to 


/('•;*)  = 


PSp  2JfKr)  ‘ 

(XR)2  nr 


(17) 
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The  corresponding  encircled  energy  of  an  aberration-free  PSF  for  a  circular  pupil  is 
given  by23 

E(rc)=  \-JQ\nrc)-J,\nrc)  (18) 

Mahajan1  and  Szapiel2  have  shown  for  the  case  of  circular  pupils  with 
rotationally  symmetric  aberrations  that  the  irradiance  distribution  of  an  aberrated  PSF 
may  be  approximated.  In  particular,  the  irradiance  distribution  of  an  aberrated  PSF 
within  the  Airy  disk  is  very  nearly  described  by  an  aberration-free  PSF  scaled  by  the 

Strehl  ratio  That  is,  within  the  Airy  disk,  the  diffraction  PSF, 

1 

/(/•)  =  4  Jexp[/cb(p)]./0(7trp)pdp  (19) 

0 

can  be  approximated  by 

r  i2 

I(r)  =  S  -  r }  ,  (20) 

7t  r 

where  S  is  the  Strehl  ratio  (given  in  Eq.  (9)  or  approximated  in  Eqs.  (11)  and  (12))  and 
the  aberration-free  PSF  has  been  normalized  to  unity  at  the  center  r  =  0  by  the  central 
irradiance  PSp/(k  R)2  The  corresponding  encircled  energy  may  similarly  be 
approximated  Mahajan  has  shown  for  the  orthonormal  Zernike  polynomials  of  fourth, 
sixth,  and  eighth  order  spherical  aberrations  (which  are  balanced)  that  the  encircled 
energy  within  small  radii  (rc  <  1  \  F)  can  be  estimated  with  a  percentage  error  of 
<  10%  for  Strehl  ratios  S  >  0. 1 .  Szapiel  provides  an  analytical  basis  for  Eq.  (20)  using 
a  Dini-sampling  method  which  involves  expressing  the  Hankie  transform  in  Eq  (19)  in 
terms  of  samples  taken  at  successive  points  given  by  the  first  few  zeros  of7,(x).2  The 
advantage  of  the  Szapiel  approach  is  that  rotationally  symmetric  aberrations  which  are 
not  balanced  tend  to  spread  the  PSF  over  a  larger  radius  and  may  be  accounted  for  by 

r  -|2 

cr(nr),  (21) 

nr 
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where  a  (n  r)  is  the  spreading  factor  (which  is  not  related  to  the  standard  deviation  of 
the  aberration).  Szapiel  proves  that  for  the  special  case  ofZemike  polynomial 
aberrations,  the  spreading  factor  <j(n  r)  =  1  The  disadvantage  of  the  Szapiel 
approach  is  that  the  calculation  of  the  spreading  factor  is  somewhat  tedious  (although 
it  can  be  done  with  a  desk-top  calculator)  and  overestimates  the  central  irradiance, 

7(0),  for  S<0.4.2 

The  focus  of  this  thesis  is  to  report  on  the  extension  of  Mahajan's  model,  Eq 

(20),  to  annular  pupils  with  rotationally  symmetric  and  nonsymmetric  aberrations,  with 

uniform  and  Gaussian  illumination  Thus,  the  general  expression  of  the  aberrated 

diffraction  PSF  given  in  Eq  (3)  is  approximated  for  spot  radii  within  the  Airy  disk  by 

/(r,e,.s)  =  S/(r,ei.e)o  =  0,  (22) 


and  for  the  uniformly  illuminated  case  by 

(23) 

nr  ner 

Similarly  the  corresponding  encircled  energies  of  aberrated  PSFs  are  also 
approximated  by  scaling  the  aberration-free  encircled  energy  by  the  Strehl  ratio 
This  thesis  reports  on  the  accuracy  of  such  a  model  for  approximating  aberrated  PSFs 
over  a  range  of  obscuration  ratios  for  primary  and  balanced  primary  aberrations 


I(r,Q,.e)  =  S ■ 


PS.. 


( XR):  (1  -e:): 
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CHAPTER  3 


NUMERICAL  RESULTS 

General  Remarks  and  Methodology 
Defocus,  spherical,  astigmatism,  coma,  and  the  corresponding  balanced 
aberrations  with  uniform  and  Gaussian  amplitudes  are  considered  in  this  thesis  The 
Gaussian  amplitude  is  given  by  A(p )  =  ^0exp(-yp2)  where  /f,,is  a  constant  and  the 
truncation  is  defined  such  that  y  =  (anf.  As  in  Eq.  (3),  a  is  the  radius  of  the  exit 
pupil,  and  w  is  the  radial  distance  at  which  the  amplitude  falls  off  to  1/e  of  the  value  at 
the  center  (known  as  the  beam  waist  of  Gaussian  laser  beams).  A  Gaussian  amplitude 
with  a  truncation  factor  of  y  -  1  is  considered. 

The  methodology  consisted  of  varying  the  obscuration  raiio,  s,  from  0,  0  1, 
0.2,...,  to  0.9.  The  value  for  o+(/l;,e)  was  estimated  using  the  Gaussian  approximation 
of  the  Strehl  ratio,  Sz  =  exp[-(o^)2],  to  obtain  Strehl  ratios  of  0  8,  0  6,  0  4,  0  2,  and 
0  1  The  magnitude  of  was  then  iteratively  refined  over  the  Gaussian 
approximation,  such  that  the  desired  Strehl  ratios  were  obtained  The  aberration 
coefficient,  Ah  was  determined  through  the  standard  deviation,  a^/l^e),  and  used  to 
calculate  the  PSF  and  encircled  energy  (which  were  normalized  by  the  central 
irradiance  PSr/(XJi )2  of  the  aberration-free  PSF)  for  each  case  of  A,  and  e  Each 
aberrated  PSF  was  then  normalized  to  unity  at  the  center  by  a  scale  factor  equal  to  S, 
i  e.,  the  PSF  was  divided  by  its  respective  Strehl  ratio,  which  by  definition  is  the  ratio 
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of  the  central  irradiances  of  the  PSF,  as  shown  in  Eq.  (9).  Encircled  energies  were 
again  calculated  for  the  PSFs  normalized  by  the  Strehl  ratio. 

The  numerical  results  presented  in  this  thesis  describe  the  difference  between 
actual  aberrated  encircled  energies  Ea  and  unaberrated  encircled  energy  Eu  scaled  by 
the  Strehl  ratio  S  The  difference  is  quantified  in  terms  of  a  percentage  error  defined 
as 


%ErrorXEXC  =  100  1 


(SXEU) 


(24) 


where  the  subscript  NENC  specifies  normalized  encircled  energy  (by  the  Strehl  ratio) 
Instead  of  using  the  actual  Strehl  ratio  of  an  aberrated  PSF  for  normalization  of 
the  aberration-free  encircled  energy,  the  Gaussian  approximation  of  the  Strehl  ratio,  Sv 
may  also  be  used  (in  practice,  the  aberratio.  variance  may  be  more  readily  measured, 
and  thus  Strehl  ratio  estimated  with  Sg ,  than  a  measurement  of  the  actual  Strehl  ratio) 
The  %oErrorSESC  in  Eq.  (24)  above  should  then  be  adjusted  for  the  percentage  error  of 
the  Strehl  approximation, 

%Errors=  ]  00(  1  -Sg  S),  (25) 

such  that 

%ErrorCaussian  =  ( %ErrorXEXC )  +  ( %ErrorSf, )  -  ( %ErrorXEXC)( %ErrorSg)  (26) 

Note  that  if  Sg  >  S  then  %ErrorSg  <  0,  and  \f\%ErrorSg\  >  %ErrorXFSV  >  0  then  the 
overall  error  of  the  model  is  improved  This  factor  may  be  significant  for  S  <0.5 
where  the  Gaussian  approximation  does  overestimate  the  actual  Strehl  ratio  (coma 
being  the  only  exception,  for  s  <  0.5)  by  more  than  a  few  percent.18 

In  contrast  to  the  scaling  of  aberration-free  encircled  energies  by  the  Strehl 
ratio,  an  alternate  method  for  estimating  encircled  energies  for  rotationally  non- 
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symmetric  aberrations  that  provides  relatively  good  accuracy  at  larger  spot  radii  and 
over  a  broader  range  cf  Strehl  ratios  is  the  comparison  of  encircled  energies  of  an 
aberrated  PSF  to  aberration-free  PSF  without  any  normalization  by  the  Strehl  ratio. 
Although  encircled  energies  vary  considerably  within  the  first  Airy  disk,  the  encircled 
energies  converge  rapidly  at  larger  radii.  As  Mahajan  has  reported,  to  a  first  order,  the 
encircled  energy  for  large  radii  is  independent  of  the  aberration  and  exhibit  an 
asymptotic  behavior  :4  The  difference  between  the  actual  encircled  energies  of 
aberrated  and  aberration-free  PSFs  is  also  quantified  in  terms  of  a  percentage  error 
defined  as 


%  Error,.  V(.  =  100 


E. 


a  J 


(27) 


The  subscript  EN(  ’  refers  to  encircle  J  energy  without  normalization  by  the  Strehl 
ratio 


Just  as  encircled  energy  increases  with  radius  (up  to  the  total  energy), 
percentage  error  of  estimated  encircled  energy  relative  to  actual  encircled  energy  for 
an  aberrated  system  with  a  "iven  obscuration  and  Strehl  ratio  changes  with  radius 
For  small  radii,  a  better  estimate  of  aberrated  encircled  energy  is  given  by  the 
aberration-free  encircled  energy  scaled  by  the  Strehl  ratio  At  larger  radii,  a  better 
estimate  of  aberrated  encircled  energy  is  obtained  by  using  the  aberration-free 
encircled  energy  without  any  scaling  by  the  Strehl  ratio. 

The  results  obtained  for  the  cases  studied  indicate  that  rotationally  symmetric 
aberrations  and  balanced  coma  have  lower  %  Error NhSC  than  %  Error Esv  over  a  larger 
range  of  radii  than  rotationally  nonsymmetric  aberrations  The  radius  for  optimal 
estimation  of  encircled  energies  without  scaling  by  the  Strehl  ratio  ranges  from  rc  ^ 
0.8X  /-'(and  larger)  for  coma,  to  rc  =  1  l?v  /-'(and  larger)  for  astigmatism,  and  rc  = 
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1  3X,  F  (and  larger)  for  rotationally  symmetric  aberrations  with  no  obscuration  (e  =  0). 
For  rotationally  symmetric  aberrations  at  large  obscuration  ratios,  e  >  0.7,  the  radius 
increases  to  greater  than  rt.  ~  2  SX  F  (hereafter  radius  will  be  specified  without  the 
units  of  XF). 

In  terms  of  percentage  errors,  %  Error SEl4C  is  less  than  %  Error Esc  for 
%ErrorSKSC  <  1 1%,  <  25%,  <  43%,  <  67%,  and  <  81%  for  Strehl  ratios  of  .V  =  0.8, 

0  6,  0  .,  0.2,  and  0  1  respectively  (see  Appendix  for  discussion  on  these  percentage 
error  values)  Practically,  percentage  errors  greater  than  20%  to  30%  may  not  be  of 
interest  Similarly,  encircled  energy  may  be  of  interest  only  within  the  central  Airy  disk 
except  where  off-axis  energy  rejection  is  a  concern 

In  addition  to  reduction  of  the  total  encircled  energy  in  a  beam  by  a  factor 
equal  to  the  fraction  of  obscured  area  of  the  pupil,  1  -  e2 ,  the  increase  of  obscuration 
ratio  affects  the  PSF  and  encircled  energy  The  radius  of  the  central  bright  spot  of  the 
aberration-free  PSF  is  slightly  reduced;  for  example,  from  r  =  1  22  at  e  =  0  to  r  =  0.8 1 
at  e  =  0.9,  as  shown  in  Figure  1 .  The  other  effect  is  seen  in  the  limit  as  the  obscuration 
ratio  approaches  unity  where  the  PSF  and  encircled  energies  of  the  rotationally 
symmetric  aberrated  beam  become  nearly  identical  to  the  aberration-free  case  when 
scaled  by  the  Strehl  ratio  within  the  second  and  even  the  third  dark  ring  The 
normalized  PSF  plots  for  the  cases  of  e  =  0  and  0.9  are  shown  in  Figures  2  and  3  for 
defocus  and  spherical  aberrations  respectively.  This  is  consistent  with  the  effect  of 
obscuration  ratio  on  the  Strehl  ratio  Mahajan  has  shown  that  Strehl  ratio  increases 
with  obscuration  for  rotationally  symmetric  aberrations  and  balanced  coma,  but 
decreases  for  astigmatism,  balanced  astigmatism,  and  coma  13 

The  results  given  here  nre  a  subset  of  all  the  data  obtained  as  part  of  this  thesis 
Results  for  each  aberration  considered  are  given  below  in  tabular  form  for  e  =  0,  0.3 

14 


I(r)/S  and  E(rcyS  l(r)  and  E(rc) 


I(r)  and  E(ij) 


l(r)/S  and  E(rcyS  l(r)  and  E(rc) 


Spheric*!,  c  •»  0 


PSFs  and  encircled  energies 
Figure  3a. 


Spherical,  t  -  0 


Normalized  PSFs  and  encircled  energies 

Figure  3b. 
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PSFs  and  encircled  energies 
Figure  3c. 


0  0  0  5  1  0  1  5  2  0  2  5  3  0  3  5  4  0 

f,  h 

Normalized  PSFs  and  encircled  energies 

Figure  3d. 
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0  5,  and  0  7.  Plots  of  PSFs,  normalized  PSFs,  encircled  energies,  and  percentage  error 
for  e  =  0  and  0,5  are  also  included  in  the  appendix  for  aberrations  other  than  defocus. 

Defocus 

The  amount  of  defocus  Ad  was  determined  by  using  the  corresponding 
standard  deviation  for  defocus  given  in  Table  1 .  Strehl  ratios  of  0.8,  0.6,  0.4,  0.2,  and 
0  1  were  obtained  for  standard  deviations  of  crM,  =  0.0744X,  0.1 1 IX,  0. 145X,  0.1 86X, 
and  0.21 3X  respectively.  The  Strehl  ratios  varied  by  less  than  1%  as  the  obscuration 
was  changed  from  e  =  0  to  0.9  for  each  corresponding  value  of  standard  deviation 
That  is,  the  Strehl  ratio  was  held  relatively  constant  about  the  desired  values  such  that 
it  was  not  necessary  to  change  the  corresponding  values  of  oM  This  is  consistent  with 
the  Gaussian  approximation,  Eq  (12),  where  Strehl  ratio  is  a  function  of  only  the 
aberration  variance.  This  implies  that  Aa r  decreases  as  e  increases  for  e  -  0  to  0.9. 

Effectively,  increasing  the  obscuration  increases  the  optical  tolerance.  The  increase  in 
depth  of  focus  by  increasing  the  obscuration  ratio  is  well  known.25 

The  percentage  difference  between  actual  encircled  energy  for  an  aberrated 
PSF  and  the  scaled  encircled  energy  of  an  aberration-free  PSF  for  circular  pupils  is  as 
low  as  <  5%  fori’>  0.1  and  rc  <  O.5.1  For  annular  pupils,  as  the  obscuration  ratio 
increases,  the  model  provides  even  better  agreement  as  shown  in  Figures  4  and  5.  For 
e  =  0.5,  S>  0. 1  and  rc  <  0.5  the  difference  is  <  4%.  Fore  =  0.7,  S  >  0. 1  and  rc  <  0.5 
the  difference  is  <  2%.  The  percentage  error  for  S  >  0  4  and  rc  <  0  8  ranges  from  < 

8%  for  e  =  0,  <  5%  for  e  =  0.5  to  <  3  %  for  e  =  0.7.  The  limiting  case  of  e  =  0.9  has  a 
difference  of  <  4%  for  S  >  0.4  and  <  14%  for  S  >  0  1  at  rc  =  2.9,  the  3rd  dark  ring. 

The  aberrated  PSF  thus  becomes  increasingly  similar  to  the  aberration-free  PSF  when 
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scaled  by  the  Strehl  ratio  for  e  =  0,9  such  that  even  for  S  =  0.1,  the  percentage  error  is 
<  2%  for  rc<  0.8 1 .  Plots  of  PSFs,  encircled  energies,  and  percentage  errors  are 
shown  in  Figures  2,  4,  and  5  for  uniformly  illuminated  pupils  and  in  Figures  6  and  7  for 
the  Gaussian  illumination  case  Percentage  errors  {%  Error jVEVC)  for  a  given  Strehl 
ratio  may  be  directly  read  off  the  plots  for  a  specified  radius  (Figures  4,  5c,  6c,  and 
7c).  Table  2  summarizes  the  percentage  errors  for  e  =  0,  0.3,  0.5,  and  0.7. 

The  Gaussian  approximation  for  the  Strehl  ratio  behaves  well  throughout  all 
values  of  e  for  the  case  of  defocus.  Percentage  error  is  relatively  unaffected  by 
variation  of  e  with  errors  ranging  from  <  9%  for  S  >  0.4  to  <  29%  for  S  >  0.2. 
Percentage  error  here  is  defined  in  Eq.  (25).  The  percentage  error  varied  by  less  than 
2%  for  a  given  Strehl  ratio  as  the  obscuration  changed  from  e  =  0  to  e  =  0  9 

Spherical  Aberration 

As  in  the  case  of  defocus,  although  the  phase  aberration  and  standard  deviation 
are  functions  of  A,  and  e,  as  can  be  seen  in  Table  1,  the  Strehl  ratio  is  largely 
dependent  on  aH.  and  not  explicitly  on  A ,  or  e  For  example,  in  the  case  of  spherical 
aberration  for  e  =  0  and  As  =  0 ,25\,  ow  -  0.07454).  and  S  =  0  8,  but  for  e  =  0.9  and  As 
=  0.25),  oM,  =  0.05885).  and  S  =  0,976.  Conversely,  for  e  =  0.9  and  ow  -  0.0754),  As 
=  0.7505).  and  S  =  0.8.  Thus  the  aberration  tolerance  increases  by  a  factor  of  three  for 
e  =  0.9  compared  to  that  for  e  =  0 

Plots  of  PSFs,  encircled  energies,  normalized  PSFs,  normalized  encircled 
energies,  and  percentage  errors  (as  defined  in  Eqs  (24)  and  (27)  above)  are  shown  in 
Figures  8  and  9  (pages  55-58)  for  spherical  aberration  with  a  uniformly  illuminated 
pupil.  The  percentage  errors  for  the  encircled  energy  approximation  with  scaling  by 
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Table  2.  Percentage  Errors.  %  Error vfcvc.  of  Estimated  Encircled  Energies  Relative 
to  Actual  Encircled  Energies  for  Defocus  With  Uniform  and  Gaussian  Illumination 
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the  Strehl  ratio  for  spherical  aberration,  at  a  radius  rc  <  0.8,  and  with  no  obscuration,  e 
=  0,  are  <  13%  for  S  >  0.4  and  <  25%  for  S  >  0.2.  Within  this  same  radius,  but  for  e  = 
0.5,  the  errors  are  <  9%  for  S  >  0  4  and  <  21%  for  S  >  0.2.  The  percentage  errors 
continue  to  decrease  as  the  obscuration  increases  and  are  <  4%  for  S  >  0.4  and  <  9% 
for  S  >  0.2  at  e  =  0.7. 

The  error  within  a  radius  defined  by  the  Airy  disk  (  r  -  1 .22)  is  <  38%  for  S  > 
0.4  and  is  <10%  for  S  >  0.8  at  s  =  0.  The  difference  is  <21%  for  S  >  0.4  and  is  <  5% 
for  S  >  0.8  at  e  =  0.5  (for  rc  <  1.22).  Perhaps  a  more  useful  radius  to  compare  the 
results  for  the  obscured,  e  =  0,5  case  is  the  slightly  reduced  radius  of  the  first  dark 
ring,  /•  =  1 .00.  The  percentage  errors  for  rc  <  1 .00  are  <  16%  for  S  >  0.4  and  <  3% 
for  S  >  0.8  for  e  =  0.5. 

Gaussian  apodization  does  not  have  a  significantly  large  effect  on  the 
percentage  error  results  and  tends  to  be  within  a  few  percent  of  the  results  for  uniform 
illumination.  The  results  are  shown  in  Table  3  and  in  Figures  10  a.,u  1 1  (pages  59-62). 

Balanced  Spherical  Aberration 

Similar  results  are  obtained  for  the  case  of  balanced  spherical  aberration.  It 
can  be  seen  from  Table  1,  that  the  balancing  defocus  for  spherical  aberration  is  given 
by  -■■■-  -( 1  +e2)^J  .  This  amount  of  defocus  minimizes  the  aberration  variance  and 
hence  results  in  the  maximum  Strehl  ratio  for  small  aberrations.26 

The  values  for  e  =  0  for  a  given  Strehl  ratio  are  the  same  values  as  used 
for  e  =  0  9  Furthermore,  the  Strehl  ratios  are  essentially  constant  as  e  is  changed  for 
each  given  value.  This  implies  that  --  exp[-(27tc^)2]  approximates  the  Strehl  ratio 
with  similar  results  for  both  circular  and  annular  pupils.  The  percentage  error  for 
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Table  3.  Percentage  Errors.  %  Error SFNC  ,_of  Estimated  Encircled  Energies  Relative 
to  Actual  Encircled  Energies  for  Spherical  Aberration  With  Uniform  and  Gaussian 
Illumination 


Uniform  Illumination  Gaussian  Illumination 
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relative  the  actual  Strehl  ranges  from  4%  to  5.2%  for  S  =  0.4,  for  an  obscuration 
0  $  e  <  0.9.The  maximum  agreement  between  the  aberrated  PSF  and  the  aberration- 
free  PSF  scaled  by  the  Strehl  ratio  is  obtained  for  the  case  of  balanced  spherical 
aberration.  As  reported  by  Mahajan,  for  e  =  0,  the  percentage  error  is  <  10%  for  S  > 

0  1  and  rc  <  1.0.  *  The  difference  increases  only  to  <  12%  at  rc<  1 .4.  The  percentage 
error  fluctuates,  it  is  negative  for  rc  <  13  then  increases  rapidly  as  radius  increases. 

The  difference  is  <  6%  for  S  >  0.2  and  <  3%  for  S  >  0.4  at  rc<  1 .3. 

As  in  the  case  of  spherical  aberration  and  defocus,  the  model  gives  more  and 
more  accurate  results  as  the  obscuration  increases.  For  e  =  0.3,  the  difference  is  <  8% 
for  5^0.1  and  <  4%  for  S'  >  0.2  at  rc<  1.3.  Fore  =  0.5  and  e  =  07  the  difference  is  < 
5%  and  <  3%  respectively  for  S  >  0. 1  at  rc  <  1 .2.  At  larger  radii,  for  e  =  0  3  at  rc  < 

2  4  the  difference  is  <31%  for  .S' >  0.4,  and  for  e  =  0.5  the  difference  is  <  13%  for  S> 
0  4  at  rc  <  2  0,  which  are  the  respective  radii  of  the  second  dark  ring  For  e  =  0.3  the 
percentage  difference  within  the  third  dark  ring  at  r  =  3.22  is  <  20%  for  S  >  0.2  and  < 
10%  f0r  S  >  0  4  The  percentage  errors  are  summarized  in  Table  4  Note  the  error 
shovvn  is  the  maximum  within  the  specified  radius.  The  plots  of  normalized  encircled 
energy  share  a  common  value  at  rc=  1.2  such  that  the  percentage  error  at  this  radius  is 
zero  as  is  shown  in  Figures  12  and  13  (pages  63-70).  Optimum  balancing  of  spherical 
aberration  with  Gaussian  apodization  is  not  in  general  given  by  the  same  amount  of 
defocus  as  for  uniformly  illuminated  pupils  Rather,  the  minimum  of  the  aberration 
variance  is  modified  to  include  the  Gaussian  weighting.9  The  results  with  Gaussian 
apodization  with  y  =  1  are  very  similar  to  uniform  illumination  and  are  shown  in 
Figures  l2e-gand  13e-g  (pages  65-66  and  69-70). 
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Table  4  Percentage  Errors.  %  Error s^,c .  of  Estimated  Encircled  Energies  Relative 
to  Actual  Encircled  Energies  for  Balanced  Spherical  Aberration  With  Uniform  and 
Gaussian  Illumination 


Uniform  Illumination  Gaussian  Illumination 
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Rotationallv  Nonsvmmetric  Aberrations 


In  general,  an  aberrated  PSF  for  astigmatism,  balanced  astigmatism,  and  coma 
when  scaled  by  S'"1  was  substantially  less  in  agreement  with  the  aberration-free  PSF 
than  the  rotationally  symmetric  aberrations.  Furthermore,  the  agreement  of  the  model 
is  largely  unimproved  as  the  obscuration  increases.  This  is  consistent  with  the  effect  of 
a  change  in  obscuration  ratio  on  the  Strehl  ratio  for  astigmatism,  balanced  astigmatism, 
and  coma  which  has  been  reported.15  The  Strehl  ratio  for  a  given  amount  of 
aberration  decreases  as  obscuration  increases  for  these  three  aberrations.  Note  that 
the  PSFs  and  normalized  PSFs  considered  here  are  plotted  along  the  Q,  =  0  axis  for 
which  the  effects  of  small  astigmatism  and  coma  aberrations  are  most  apparent 
(assuming  cos0  terms  in  the  aberration  and  no  sin0  terms). 

Astigmatism 

The  percentage  difference  at  rc  <  0.8  ar.d  e  =  0  is  <  5%  for  S  >  0.8,  <  1 1%  for 
S  >  0.6,  and  <  20%  for  S  >  0.4.  The  error  does  not  change  by  more  than  2%  for 
larger  obscurations  up  to  e  =  0.7  within  the  radius  of  rc  <  0.8.  The  error  more  than 
doubles  at  rc  <  1 .2  and  is  <  11%  for  S  >  0.8,  <  25%  for  S  >  0.6,  and  <  52%  for  S  >  0.4 
at  e  =  0  and  fluctuates  by  5%  as  the  obscuration  increases  (S  =  0.8)  and  up  to  12%  (S 
=  0.4)  for  obscurations  e  <  0.7.  Thus,  although  the  percentage  error  does  not  change 
considerably  with  a  change  in  obscuration  ratio,  the  model  has  a  limited  applicability 
to  large  Strehl  ratios  of  about  S  >  0.4  for  errors  in  the  order  of  20%  or  less  at  a  radius 
rc  <  0  8  Plots  of  PSFs,  normalized  PSFs,  encircled  energies,  and  percentage  error  are 
shown  in  Figures  14  and  15  (pages  71-74)  for  uniform  illumination  and  Figures  16  and 
17  (pages  75-78)  for  Gaussian  illumination. 
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Alternatively,  the  percentage  differences  between  the  actual  encircled  energies 
of  aberrated  and  aberration-free  PSFs  without  normalization  by  the  Strehl  ratio 
provide  a  better  estimate  for  radii  rc>  1.2  For  e  =  0  at  rc  =  1.2,  the  percentage 
differences  (now  defined  in  Eq.  (27)  and  are  negative)  between  actual  aberrated  and 
aberration-free  encircled  energies  are  10%,  25%,  50%,  and  450%  for  S  =  0  8,  0  6,  0.4, 
and  0  2,  respectively.  The  error  drops  to  3%,  5%,  14%,  and  37%  for  the  same 
respective  Strehl  ratios  at  rc  =  2.0.  At  rc  >  3.0,  the  percentage  difference  for  S  >  0  2  is 
<  9%.  Table  5  summarizes  the  results 

Since  one  of  the  effects  of  an  annular  pupil  is  to  increase  the  energy  in  the 
outer  regions  relative  to  the  energy  within  the  Airy  disk,  it  is  to  be  expected  that  for 
larger  obscurations,  the  encircled  energies  do  not  converge  as  rapidly.  For  e  =  0.3  at 
rc  =  2.0  the  percentage  differences  are  4%,  12%,  22%,  and  42%  for  S  =  0.8,  0.6,  0.4, 
and  0  2  respectively  The  encircled  energies  at  rc  >  3.0  converge  to  the  same 
percentage  differences  as  for  e  =  0,  <  5%.  The  results  for  e  -  0.5  are  within  2%  of  the 
values  for  e  =  0  3  at  rc  >  2.0.  The  most  significant  difference  for  e  =  0.5  is  that  the  S  = 
0  1  case  is  within  27%  of  the  aberration-free  encircled  energy  for  radius  of  rc  =  3.0. 
The  percentage  differences  for  e  >  0.7  at  rc  =  2  0  are  3%,  8%,  16%,  and  28%  for  S  -■ 

0  8,  0  6,  0.4,  and  0.2  respectively.  For  the  highly  obscured  case  of  e  =  0  9,  the 
percentage  difference  is  <  1 0%  for  S  >  0. 1  at  a  radius  of  rc  >  2  1 

Although  encircled  energies  at  regions  greater  than  the  Airy  disk  may  not  be  of 
interest  in  some  applications,  some  useful  generalizations  can  be  made  in  the  case  of 
astigmatism.  The  region  where  a  better  estimate  of  encircled  energy  is  obtained  by 
using  actual  aberration-free  encircled  energy  rather  than  that  scaled  by  the  Strehl  ratio 
is  a  radius  of  rc  >  1 .2  for  S  >  0.6.  The  same  is  true  for  a  radius  of  rc  >  1.3  for  S  >  0.4 
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Table  5.  Percentage  Errors.  %  Error >,£>r.  of  Estimated  Encircled  Energies  Relative 
to  Actual  Encircled  Energies  for  Astigmatism  With  Uniform  and  Gaussian  Illumination 
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and  rc  >  1 .5  for  S  >  0.2.  These  regions  are  valid  for  all  cases  of  obscuration  ratios  e  = 
0  to  e  =  0.9  as  shown  in  Figures  14d  and  15d  (pages  72  and  74). 

Balanced  Astigmatism 

The  results  obtained  for  balanced  astigmatism  are  similar  to  astigmatism 
Here,  balanced  astigmatism  is  defined  as  astigmatism  with  a  balancing  defocus,  where 
Ad  =  -t  Aa  is  introduced  into  the  aberration.  The  percentage  errors  given  in  Table  6 

refer  to  the  percentage  difference  between  actual  encircled  energies  of  aberrated  PSF 
and  aberration-free  encircled  energy  scaled  by  the  Strehl  ratio.  For  e  =  0,  the  errors 
within  rc  <  0.8  are  <  6%,  <13%,  and  <  23%  for  S  >  0.8,  0.6,  and  0.4,  respectively.  At 
rc<  1.2  the  errors  are  <  14%  for  S  >  0  8  and  <  28%  for  S  >  0  6  The  percentage 
errors  for  e  =  0.3  are  the  same  as  at  e  =  0  with  the  exception  that  for  S  >  0.6  and  S  > 
0.4  the  error  is  <  12%  and  <  22%  respectively  at  rc  <  0.8.  Similarly,  the  percentage 
errors  are  within  2%  of  values  at  e  =  0  for  e  =  0.5.  For  e  =  0.7  the  percentage  errors 
are  different  than  e  =  0  only  for  rc<  1.2  and  are  <  16%  and  <  34%  for  5  >  0.8  and  S  > 
0.6  respectively  The  results  are  summarized  in  Figures  18  through  21  (pages  79-86) 
and  Table  6  for  uniform  illumination  and  Gaussian  illumination  with  a  truncation  factor 
of  y  =  1 

Following  the  same  approach  as  astigmatism  above,  the  percentage  difference 
of  encircled  energy  for  an  aberrated  PSF  relative  to  an  aberration-free  PSF  without 
scaling  by  the  Strehl  ratio  is  significantly  less  than  that  scaled  by  the  Strehl  ratio  at 
larger  radii  The  percentage  differences  of  actual  encircled  energies  relative  to 
aberration-free  encircled  energy  for  e  =  0,  are  <  2%,  <  4%,  <  7%,  and  <  13%  for  S  > 

0  8,  0.6,  0.4,  and  0.2  respectively  at  a  radius  rc  >  2.0.  For  e  =  0.3  and  e  -  0  5  at  rc  > 
2,0  the  differences  are  <  3%,  <  7%,  <  13%,  and  <  24%  for  the  same  values  of  Strehl 
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Table  6  Percentage  Errors,  %  Error SESC.  of  Estimated  Encircled  Energies  Relative 
to  Actual  Encircled  Energies  for  Balanced  Astigmatism  With  Uniform  and  Gaussian 
illumination 


Balanced  Astigmatism  Balanced  Astigmatism 

Uniform  Illumination  Gaussian  Illumination 


s  > 

rc<  0  8 

rc<  1.0 

rc<  1.2 

e  =  0 

rr<  0.8 

rc<  10 

rc'- 

08 

6 

10 

14 

6 

9 

12 

0.6 

13 

20 

28 

14 

20 

26 

0.4 

23 

34 

44 

23 

33 

41 

0  2 

38 

52 

64 

e  -  0.3 

36 

50 

60 

0.8 

5 

9 

14 

5 

9 

14 

0.6 

12 

20 

27 

12 

J  9 

26 

0.4 

22 

33 

44 

21 

31 

41 

0.2 

36 

51 

64 

e  =  0.5 

35 

49 

60 

0.8 

5 

9 

14 

5 

9 

14 

0.6 

10 

20 

30 

10 

19 

28 

0.4 

20 

34 

46 

19 

32 

44 

0.2 

35 

54 

66 

e  =  0.7 

33 

51 

65 

0.8 

5 

11 

16 

5 

10 

16 

0  6 

1 1 

24 

34 

11 

29 

34 

0.4 

nn 

40 

52 

21 

39 

51 

0.2 

40 

61 

73 

39 

60 

71 

ratio.  For  e  =  0.7  the  percentage  differences  are  slightly  improved:  <  3%,  <  6%, 

<  1 1%,  and  <  19%.  At  a  smaller  radius  of  rc>  1.2  the  percentage  differences  are 

<  8%  and  '  20  %  for  5  >  0.8  and  S  >  0.6  respectively,  for  all  obscurations  0  <  e  <  0.9, 
which  is  a  smaller  error  than  that  obtained  with  the  scaled  encircled  energies 
considered  above 

Coma 

A  coma  aberration  is  known  to  cause  a  line  of  sight  error  which  is  essentially  a 
shift  of  the  centroid  of  the  PSF  away  from  the  central  axis 26  AJthough  Mahajan  has 
shown  that  the  peak  and  centroid  of  the  PSF  are  not  coincident  for  an  aberrated 
system,  the  peak  also  shifts  with  a  coma  aberration,  p"cos0,  as  well  as  any  Zemike 
aberration  with  7^(p)cos0  terms  were  n  is  an  odd  integer.26  Given  that  the  peak  of 
the  PSF  is  displaced  from  the  origin,  the  Strehl  ratio  does  not  describe  the  maximum 
value  of  the  PSF,  but  rather  the  irradiance  at  the  center.  For  example,  for  an 
aberration  of  Ac  =  0  8 1 X  and  an  obscuration  ratio  of  e  =  0,  the  Strehl  ratio  is  5  =  0 1 0 
whereas  the  maximum  value  of  the  PSF  is  0.7  (normalized  to  the  e  =  0,  aberration-free 
PSF)  and  is  located  a  distance  r  =  1 .05  away  from  the  central  axis.  For  an  obscuration 
of  e  =  0  9  and  an  aberration  of  Ac  =  0.34)w,  the  Strehl  ratio  is  still  5  =  0.10  although 
the  PSF  maximum  is  0  99  (normalized  to  the  e  =  0  9,  aberration-free  PSF)  and  is 
located  at  a  distance  r  -  0.60  away  from  the  center.  It  has  been  reported  that  the 
distance  of  the  peak  from  the  origin  increases  monotonicallv  with  Ac  up  to  Ac  ^  1  6X 
for  e  =  0  and  up  to  Ac  =  2.5X  for  e2  =  0.5,  but  fluctuates  for  larger  aberrations.26 

For  a  given  value  of  aberration  variance,  as  the  obscuration  ratio  increases,  the 
Strehl  ratio  decreases  Conversely,  the  peak  of  an  aberrated  PSF  increases,  relative  to 


Table  7.  Percentage  Errors  for  Coma.  With  and  Without  Strehl  Ratio  Scaling,  for 
Uniform  Illumination  Note  That  Results  With  Gaussian  Apodization  Are  Within 
2%  of  the  Values  Given  Below 


%ErrorNKSC 

%ErrorESV 

Si 

rc<  0.8 

rc<  '  2 

rc>\  2 

rc>  2.0 

rc-  2-5 

e  =  0 

0  8 

11 

18 

4 

3 

1 

0.6 

24 

36 

8 

5 

2 

0.4 

37 

53 

17 

10 

3 

e  =  0.3 

0.8 

12 

18 

5 

5 

1 

0.6 

25 

37 

12 

10 

2 

0.4 

38 

54 

20 

18 

4 

e  =  0.5 

0.8 

13 

21 

6 

4 

2 

0.6 

27 

40 

13 

9 

3 

0.4 

43 

58 

23 

22 

5 

II 

© 

'-4 

0.8 

16 

23 

5 

3 

3 

0.6 

33 

43 

12 

7 

7 

0.4 

50 

62 

8 

10 

10 

an  aberration-free  PSF,  with  the  obscuration  ratio  for  a  given  value  of  aberration 
variance.  Additionally,  the  widfh  of  the  PSF,  or  the  "spot  size,"  was  relatively 
unchanged  for  a  given  obscuration  ratio  for  the  aberration  considered  (S  >  0. 1)  as 
shown  in  Figures  22a  and  23a  (pages  87  and  89)  This  suggests  that  with  an  optimal 
amount  of  tilt  introduced,  the  model  presented  in  this  thesis  provides  high  accuracy  for 
balanced  coma  aberration.  Just  as  in  the  case  of  astigmatism  and  balanced  astigmatism 
above,  the  percentage  error  of  the  encircled  energy  of  the  aberration-free  beam  scaled 
by  the  Strehl  ratio  relative  to  the  actual  encircled  energies  of  a  beam  aberrated  with 
coma  are  considerably  larger  than  the  results  obtained  for  rotationally  symmetric 
aberrations  For  S  >  0  8  at  rc  <  0  8  the  error  is  <  1 1%  for  e  =  C  and  increases  steadily 
to  <  16%  for  e  =  0.7.  Similarly,  for  S  >  0.6  at  rc  <  0.8  the  error  is  <  24%  for  e  =  0 
and  increases  to  <  33%  at  e  =  0.7.  The  piots  for  PSFs,  normalized  PSFs,  encircled 
energies,  and  percentage  error  are  shown  in  Figures  22  through  25  (pages  87-94) 

The  percentage  errors  are  summarized  in  Table  7. 


as  can  be  seen  from  the  phase  aberration  in  Table  1 .  As  with  balanced  spherical 
aberration,  A,  is  determined  such  that  aberration  variance  is  a  minimum  (which  is 
affected  by  Gaussian  apodization  as  before)  This  amount  of  tilt  maximizes  the  Strehl 
ratio  for  small  aberrations,  Ac  <  0.7/,. 26  The  tilt,  given  by  P  =  At  a  where  a  is  the 
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outer  radius  of  the  pupil  in  units  of  waves,  is  the  amount  of  angular  or  pointing 
compensation  for  the  line  of  sight  error  due  to  coma. 

In  general,  for  the  cases  studied,  the  effect  of  introducing  a  tilt  given  by  Eq 
(28)  was  consistent  in  over-correcting  the  effect  of  coma  and  shifting  the  peak  beyond 
the  origin  (in  the  opposite  direction  as  the  original  coma)  for  larger  aberrations.  This 
can  be  seen  in  Figures  26b  and  27b  (pages  95  and  97).  The  maximum  value  of  the 
PSF  is  centered  at  the  origin  (within  the  sampling  resolution  where  the  sample  size  is 
0.05A  F)  for  Strehl  ratios  S  >0,8  which  corresponds  to  Ac  =  0.64^.  It  is  evident  that 
at  Strehl  ratios  less  than  S  -  0.8  the  amount  of  tilt  given  by  Eq.  (28)  does  not 
maximize  the  Strehl  ratio  for  e  =  0  At  larger  obscuration  ratios,  the  tilt  given  in  Eq 
(28)  more  effectively  balances  the  coma  by  moving  the  peak  of  the  PSF  closer  to  the 
center  for  larger  amounts  of  coma  aberration.  At  e  =  0.7,  A,  given  in  Eq.  (28)  centers 
the  peak  of  the  PSF  for  S  >  0.4  which  is  for  a  coma  o(Ac  =  1.7  IX. 

Unlike  the  rotationally  symmetric  aberrations  discussed  above,  the  percentage 
errors  for  balanced  coma  do  not  change  considerably  (<  4%  for  S  >  0.2)  as  the 
obscuration  ratio  increases  for  radii  rc  <  1 .2.  For  larger  radii,  however,  the  errors 
decrease  as  obscuration  increases  such  that  at  e  =  0.7,  the  error  is  approximately  one- 
half  that  at  e  =  0  at  radii  rc  >  1 .65  (error  for  S  =  0.4  is  21%  at  e  =  0  7).  The 
percentage  errors  obtained  are  comparable  to  rotationally  symmetric  aberrations  as 
shown  in  Figures  26  through  29  (pages  95-102).  The  percent  errors  are  given  in  Table 
8.  At  rc  <  0.8  and  e  =  0  the  errors  are  <  9%  for  S  >  0.2  and  <  1 8%  for  S  >  0. 1 .  At  rc 
<  0.8  and  e  =  0.3  and  e  =  0.5  the  errors  are  <  7%  for  S  >  0.2  and  <  17%  at  S  >  0. 1 . 
The  errors  increase  somewhat  for  e  =  0.7,  being  <  1 0%  and  <  24%  for  S  >0.2  and  S  > 
0.  l  respectively.  Results  for  the  Gaussian  illumination  case  are  approximately  double 
the  percentage  errors  as  for  the  uniform  illumination  case  at  radius  rc  <  0.8  but  are 
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Table  8  Percentage  Errors,  %  Error c.  of  Estimated  Encircled  Energies  Relative 
to  Actual  Encircled  Energies  for  Balanced  Coma  Aberration  With  Uniform  and 
Gaussian  Illumination 


Uniform  Illumination  Gaussian  Illumination 


s> 

rc.<0.8 

rc-  12 

8  =  0 

rc-  0  8 

rc-  1-2 

0  8 

1 

4 

2 

5 

0.6 

2 

9 

5 

12 

0.4 

4 

18 

9 

23 

0.2 

9 

37 

18 

42 

0.1 

18 

55 

o' 

II 

10 

28 

58 

0.8 

1 

3 

2 

4 

0.6 

1 

8 

4 

11 

0.4 

2 

18 

7 

22 

0.2 

7 

38 

15 

43 

0.1 

16 

59 

e  =  0.5 

28 

63 

0.8 

1 

3 

2 

4 

0.6 

1 

8 

3 

9 

0.4 

2 

18 

7 

20 

0.2 

7 

40 

15 

43 

0.1 

17 

60 

e  =  0.7 

28 

64 

0.8 

1 

2 

3 

3 

0.6 

1 

7 

5 

7 

0.4 

3 

16 

Q 

16 

0.2 

10 

40 

17 

38 

0.1 

24 

61 

31 

60 

43 


only  slightly  larger  than  uniform  illumination  at  a  radii  of  rc  *  1 .2. 


Mixed  Aberrations 

So  far  we  have  discussed  the  PSFs  aberrated  by  a  single  primary  aberration. 
Now  we  briefly  consider  examples  of  mixed  aberrations.  A  mixture  of  spherical, 
astigmatism,  and  coma  was  studied,  as  well  as  a  mixture  of  spherical,  astigmatism, 
coma,  and  defocus.  The  balancing  defocus  that  was  introduced  is  the  same  as  used  for 
optimally  balanced  spherical  aberration,  namely,  Ad  =  -(l+e2)/!^  which  minimizes  the 
aberration  variance  of  spherical  aberration  only.  The  same  standard  deviation  of 
aberration  was  us'.d  to  calculate  the  respective  aberrations,  thus,  the  magnitudes  of 
each  type  of  aberration  was  not  equal.  For  example,  for  a  Strehl  ratio  of  S  =  0. 1,  a 
peak  aberration  ofAs  =  1.5X,  Ad  =  -1.5X,  Ac  =  1.3X,  and  Aa  =  1.8X  was  used.  The 
PSFs,  encircled  energies,  and  peicentage  error  plots  are  shown  in  Figures  30  and  3 1 
(pages  103-106)  for  the  aberration  mixture  without  defocus  and  in  Figures  32  and  33 
(pages  107-1 10)  for  the  mixture  with  defocus.  The  percentage  errors  of  the  encircled 
energy  of  the  aberration-free  PSF  scaled  by  the  Strehl  ratio  relative  to  the  actual 
encircled  energy  of  the  aberrated  PSF  are  shown  in  Table  9.  Note  the  percentage 
errors  are  very  comparable  to  errors  for  cases  of  the  other  rotationally  nonsymmetric 
aberrations,  with  the  exception  of  balanced  coma  which  behaved  similar  to  the 
rotationally  symmetric  aberrations. 
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Table  9,  Percentage  Errors.  %  Error SEKCi_ of  Estimated  Encircled  Energies  Relative 
to  Actual  Encircled  Energies  for  Mixed  Aberrations  With  Uniform  Illumination 


Spherical,  Coma, 

Spherical,  Coma,  Astigmatism  Astigmatism,  &  Defocus 


s> 

f.  =  0 

rc<0.8 

r,<1.0 

rc<1.2 

S> 

6  =  0 

rc<0.8 

rc  <  1 .0 

rc  '■ 

0.87 

4 

5 

8 

0.8 

9 

12 

15 

0.64 

11 

17 

23 

0.6 

17 

24 

28 

0.43 

20 

29 

36 

0.4 

20 

28 

35 

0.20 

29 

41 

50 

0.2 

28 

40 

50 

0.09 

e  =  0.3 

44 

58 

67 

0.1 

e  =  0.3 

35 

50 

62 

0  87 

4 

5 

7 

0.8 

8 

12 

15 

0.64 

11 

16 

20 

0  6 

18 

25 

30 

0.42 

20 

28 

35 

0.4 

24 

32 

39 

0.19 

27 

37 

45 

0.2 

30 

43 

54 

0.08 

e  =  0.5 

40 

54 

64 

0.1 

e  =  0.5 

39 

55 

66 

087 

4 

5 

6 

0  8 

9 

14 

16 

0.65 

10 

16 

20 

0.6 

20 

28 

34 

0.42 

20 

28 

34 

0.4 

31 

42 

49 

0  17 

24 

34 

41 

0.2 

35 

49 

60 

0.07 

6  =  0.7 

31 

45 

55 

0.1 

6  =  0.7 

44 

59 

70 

088 

5 

8 

9 

0.8 

11 

16 

19 

0.66 

12 

18 

22 

0.6 

24 

34 

37 

042 

23 

33 

37 

0.4 

39 

50 

55 

0,13 

32 

43 

49 

0.2 

46 

62 

71 

005 

34 

46 

51 

0.1 

55 

67 

72 
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CHAPTER  4 


DISCUSSION  AND  CONCLUSIONS 

The  aberrated  point-spread  function  and  encircled  energy,  may  be 
approximated  within  the  central  Airy  disk  for  both  circular  and  annular  pupils  for  each 
of  the  primary  and  balanced  primary  aberrations  using  the  aberration-free  PSF 
normalized  by  the  Strehl  ratio.  The  range  of  spot  radii  that  may  be  used  in  this 
approximation  depends  on  the  amount  and  type  of  the  aberration.  The  effect  of  a 
central  obscuration  on  the  accuracy  of  the  model  also  depends  on  the  type  of 
aberration  For  defocus,  spherical,  balanced  spherical,  and  balanced  coma  aberrations, 
an  increase  in  obscuration  ratio  improves  the  accuracy  of  the  model  (i  .  e.,  the 
percentage  error  decreases)  An  increase  in  obscuration  ratio  has  little  effect  on 
accuracy  for  astigmatism  and  balanced  astigmatism  for  S  >  0.4  and  tends  to  increase 
the  percentage  error  for  smaller  Strehl  ratios  Coma  without  a  balancing  tilt  is  the  only 
aberration  for  which  the  accuracy  of  the  model  degrades  significantly  with  an  increase 
in  obscuration  ratio.  With  the  exception  of  coma  and  balanced  astigmatism,  the 
percentage  error  of  the  estimated  encircled  energy  relative  to  the  actual  value  is  no 
greater  than  20%  for  Strehl  ratios  S  >  0.4  at  spot  radii  rc  <  0.8  X  F. 

The  model  accuracy  is  maintained  for  the  mixed  aberrations  of  spherical, 
astigmatism,  and  coma  at  rc  <  0 .8  X  F,  with  slightly  better  accuracy  than  for  other 
rotationaily  nonsymmetric,  single  aberrations  (except  balanced  coma)  Mixed 
aberrations  with  a  balancing  defocus  have  larger  percentage  error  than  mixed 


46 


aberrations  without  defocus,  by  nearly  double  the  amount  Even  with  defocus, 
however,  the  percentage  errors  are  lower  than  coma  aberration  without  balancing  tilt. 

The  effect  of  a  Gaussian  apodization  (with  1/e2  irradiance  at  the  edge  of  the 
pupil  compared  to  its  central  value  in  the  absence  of  obscuration)  was  to  degrade  the 
accuracy  of  the  model  for  balanced  coma  to  approximately  twice  the  error  obtained  for 
the  case  of  uniform  illumination  for  radii  rc  <  0  8  X  F  but  was  within  5%  of  uniform 
illumination  results  at  rc  «  1.2  \F.  Gaussian  apodization  results  for  astigmatism, 
balanced  astigmatism,  coma,  and  the  rotationally  symmetric  aberrations  were  within 
5%  of  the  results  obtained  for  uniform  illumination  for  Strehl  ratios  S  >  0.4  and  for 
radii  rc<  1.0  XF. 

The  applicability  of  the  model  over  a  given  radii  is  directly  affected  by  how  the 
actual  encircled  energies  of  aberrated  optical  systems  converge  to  the  aberration-free 
encircled  energy  as  the  radius  of  the  spot  increases.  The  asymptotic  behavior  depends 
on  the  type  of  the  aberration  and  obscuration  ratio.  The  radii  for  which  a  better 
estimate  of  encircled  energy  is  obtained  (i  e  ,  less  percentage  error)  by  using  the  actual 
aberration-free  encircled  energy  rather  than  the  aberration-free  encircled  energy  scaled 
by  the  Strehl  ratio  ranges  from  rc  >  0  8  \  F  for  coma  without  tilt  to  rc  >  1.2  \  F  for 
astigmatism,  balanced  astigmatism,  balanced  coma,  defocus  (e  <  0  5),  and  spherical  (e 
<  0  5)  aberrations  For  balanced  spherical  and  other  rotationally  symmetric  aberrations 
(at  large  obscuration  ratios,  e  >  0.5),  actual  encircled  energies  converge  much  more 
slowly,  and  the  encircled  energy  is  better  estimated  with  the  aberration-free  encircled 
energy  scaled  by  the  Strehl  ratio  out  to  radii  of  at  least  the  second  and  even  third  dark 
ring. 

The  behavior  of  balanced  coma  and  rotationally  symmetric  aberrations  at  large 
obscuration  ratios  results  in  a  high  degree  of  model  accuracy  when  the  encircled 
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energy  is  estimated  using  the  aberration-free  encircled  energy  scaled  by  the  Strehl  ratio 
of  the  aberrated  PSF.  This  effect  can  be  understood  intuitively  by  analyzing  the 
behavior  of  the  aberration  variance  of  these  aberrations  in  the  limit  ofe  — >  1  (see  Table 
1)  Note  that  unlike  the  other  aberrations  (astigmatism,  balanced  astigmatism,  and 
coma)  the  aberration  variance  equals  zero  when  the  obscuration  approaches  unity. 

The  effect  of  the  aberration  then  tends  to  scale  the  central  region  of  the  PSF  by  the 
Strehl  ratio  but  otherwise  does  not  disturb  the  PSF  in  this  central  region  This 
approach  to  predicting  the  behavior  cannot,  however,  hold  true  as  the  radius  increases 
to  large  values  because  of  the  asymptotic  behavior  of  all  aberrations  wherein  the 
encircled  energy  of  even  a  highly  obscured,  aberrated  beam  must  eventually  converge 
to  the  aberration-free  encircled  energy  (i.  e.,  the  total  energy) 

Areas  where  future  work  may  be  worthwhile  include  revisiting  the  Szapiel 
approach  of  modifying  the  basic  equation  for  the  model,  Eq.  21,  to  account  for 
aberrations  that  tend  to  spread  the  spot  radii.  An  attempt  to  simplify  the  determination 
of  such  a  factor  could  permit  application  of  the  model  described  in  this  thesis  with  a 
simple  adjustment  that  would  be  used  for  certain  aberrations  (such  as  defocus  and 
spherical  aberration  at  low  obscuration  ratios  and  rotationally  nonsymmetric 
aberrations). 
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The  range  of  radii  for  optimum  estimation  approach  (that  scaled  by  Strehl 
versus  not  scaled  by  Strehl)  can  be  determined  by  setting  the  percentage  errors  defined 
in  Eqs.  (24)  and  (27)  equal  to  each  other.  Thus  for  a  given  Strehl  ratio,  S,  the  region 
where  a  better  estimate  using  the  aberration-free  encircled  energy  scaled  by  the  Strehl 
ratio  is  obtained  is  for: 

F  2 

!*-> -  (29) 

Ea  ]  +  S 

Thus  the  percentage  errors  where  %  Error SESC  is  less  than  %  Error ESC  are 
%oErrorNENC  <  1 1%,  <  25%,  <  43%,  <  67%  and  <  81%  for  Strehl  ratios  of  S  =  0.8, 

0.6,  0.4,  0.2,  and  0.1  respectively. 

In  light  of  the  asymptotic  behavior  of  encircled  energy  which  results  in  the 
encircled  energies  of  an  aberrated  system  converging  to  the  aberration-free  value  and 
eventually  to  the  total  energy  as  radius  increases,  the  asymptotic  behavior  of  the 
percentage  errors  may  similarly  be  understood.  The  following  conditions  hold  true: 

as  Ea  — >  Eu  (30) 

%  Error NESC  — ^  1 00  ( 1  -5) 

and  %  Error exc  — ^  0 

Conversely, 

as  Ea  — >  (S)(EJ  (31) 

%  Error SESC  — ^  0 
and  %  Error Esc  — ^  1 00  ( 1  - 1/5). 
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Figures  8  through  35 

Plots  of  PSFs,  normalized  PSFs,  encircled  energies,  and  percentage  errors 
for  spherical,  balanced  spherical,  astigmatism,  balanced  astigmatism,  coma, 
balanced  coma,  and  mixed  aberrations 


PSFs  and  encircled  energies 
Figure  8a. 


Normalized  PSFs  and  encircled  energies 

Figure  8b. 


PSFs  and  encircled  energies 
Figure  12a. 


Balanced  Spherical,  e  =  0 
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Figure  12b. 


Balanced  Spherical,  e  = 


I(r)/S  and  E(rc)/S  J(r)  and  E(^) 


Balanced  Spherical  with  Gaussian  Apod,  c  *  0.5 


r.  h 

PSFs  and  encircled  energies 
Figure  13e. 


Balanced  Spherical  with  Gaussian  Apod  z  =  0.5 


Normalized  PSPs  and  encircled  energies 

Figure  13f. 
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Astigmatism  with  Gaussian  apodization.  e  ■  0 
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PSFs  and  encircled  energies 
Figure  20a. 
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Normalized  PSFs  and  encircled  energie 

Figure  20b. 


Normalized  PSFs  and  encircled  energies 

Figure  21b. 
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Figure  22c. 
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Figure  22d. 


PSFs  and  encircled  energies 
Figure  23a. 


Normalized  PSFs  and  encircled  energies 

Figure  23b. 


Coma  with  Gaussian  apodization.  t  *  0 


Coma  with  Gaussian  apodization,  e  =  0, 


Coma  with  Gaussian  anodization,  e 


PSFs  and  encircled  energies 
Figure  26a. 


Normalized  PSFs  and  encircled  energies 

Figure  26b. 


PSFs  and  encircled  energies 
Figure  27a. 


Normalized  PSFs  and  encircled  energies 

Figure  27b. 
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Normalized  PSFs  and  encircled  energies 

Figure  28b. 
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PSFs  and  encircled  energies 
Figure  30a. 
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Normalized  PSFs  and  encircled  enrgies 

Figure  30b. 
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PSFs  and  encircled  energies 
Figure  31a. 
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Normalized  PSFs  and  encircled  energies 

Figure  31b. 
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PSFs  and  encircled  energies 

Figure  32a. 
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Normalized  PSFs  and  encircled  energies 

Figure  32b. 
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